function [res_coreg] = Coreg_config_auto(Y,Ks, nCluster)
[nSmp, ~, nKernel] = size(Ks);
res_coreg = [];
options = [];
options.lambda=0.01;
H = coreg_centroid_on_multi_kernel(Ks, nCluster, options);
H_normalized = H ./ repmat(sqrt(sum(H.^2, 2)), 1,nCluster);
for i1 = 1:50
clusterIdx = litekmeans(H_normalized, nCluster, 'MaxIter',100, 'Replicates', 1);
res_coreg = [res_coreg; ClusteringMeasure(Y, clusterIdx)];
end